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We apply the sum-over-states approach to calculate partial contributions to the parity non- 
conservation (PNC) in cesium [Porsev et al, Phys. Rev. D 82, 036008 (2010)]. We have found 
significant corrections to two non-dominating terms coming from the contribution of the core and 
highly excited states (n > 9, the so called tail). When these differences are taken into account the 
result of Porsev et al, Epnc = 0.8906(24) x 10""i(-Qw-/A'') changes to 0.8977(40), coming into 
good agreement with our previous calculations, 0.8980 (45). The interpretation of the PNC mea- 
surements in cesium still indicates reasonable agreement with the standard model (1.5 ct), however 
gives new constraints on physics beyond it. 

PACS numbers: ll.30.Er, 31.15.A- 



INTRODUCTION 

The search for new physics beyond standard model us- 
ing parity non-conservation (PNC) in atoms culminated 
in 1997 when Boulder group reported very accurate mea- 
surements of the PNC 6s — 7s amplitude in cesium [ij. 
The experimental uncertainty was only 0.35%. Interpre- 
tation of the measurements based on early calculations 
by Novosibirsk [3] and Notre Dame [1] groups indicated 
good agreement with the standard model. However, the 
declared theoretical accuracy of these early calculations 
(1%) didn't match the accuracy of the measurements. 
Bennett and Wieman re-analyzed the accuracy of the cal- 
culations by comparing calculated observables with new 
experimental data [J]. They pointed out that many dis- 
crepancies between theory and experiment were resolved 
in favor of theory and suggested that actual accuracy of 
the calculations was 0.4%. This lead to about 2.3 ct 
deviation of the value of the weak charge of the cesium 
nucleus from the prediction of the standard model. The 
discrepancy were resolved when Breit [^, Q and radia- 
tive (see, e.g and references therein) corrections were 
included into calculations. On the other hand, new cal- 
culations of the correlations [1, 0] didn't change the old 
results of and rather confirmed the suggestion made 
in Q that their accuracy was high. The new value of the 
weak nuclear charge was about 1 ct smaller than the pre- 
diction of the standard model which should be considered 
as good agreement. 

The situation changed when the latest, most sophisti- 
cated calculations of the PNC in Cs were reported by Por- 
sev et al [13] ■ The authors of this work used the sum-over- 
states approach and applied the coupled-cluster with sin- 
gle, double and valence triple excitations for the leading 
terms in the sum. They claimed 0.27% uncertainty of the 
calculations while their correlated PNC amplitude was 
about 0.9% smaller than in previous calculations. This 
led to perfect agreement with the standard model, with 
central points for the weak nuclear charge extracted from 



the measurements and predicted by the standard model 
coinciding exactly: QwC^^Cs) = -73.16(29)oxp(20)thcor 
and 0|^("3Cs) = -73.16(3) The smaller value of 
the correlated PNC amplitude is attributed in [10] to the 
role of higher-order correlations. 

Although all old and new calculations of the PNC in 
Cs lead to agreement with the standard model, the re- 
sults of [13] have important implications imposing strong 
constraints on new physics beyond the standard model. 
Therefore, it is worth studying further the reasons for the 
difference in the calculations. 

The authors of pXi] paid great deal of attention to the 
leading terms, performing very sophisticated calculations 
for them and demonstrating high accuracy by comparing 
with available experimental data. The uncertainty for 
the minor terms was assumed to be 10% based on the 
spread of the values in different approximations. The 
sum-over-states approach used in [10] has an important 
shortcoming. Calculation of each term in the sum is 
practically independent of others and therefore the high 
accuracy for the leading terms does not guarantee high 
accuracy for the sum. 

In this paper we use the sum-over-states approach to 
study possible reasons for the difference between the re- 
sults of [13] and previous calculations [3, Q . We assume 
that the main term was calculated correctly in [l^] and 
focus our attention on the minor terms such us contri- 
bution of the core states and liighly excited (tail) states. 
We include core polarization and correlation corrections 
into the core and tail contributions and find significant 
difference for both these terms between our calculations 
and those reported in [l^l . Our core contribution has dif- 
ferent sign while being similar in value. We have agree- 
ment with for the tail contribution when only core 
polarization effects are taken into account. However, fur- 
ther inclusion of Bruekner-type correlations increase the 
PNC amplitude beyond the 10% uncertainty for the tail 
assumed in ^1(| . 

If core and tail contributions of [13] are replaced by 
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the findings of present work the resulting PNC ampli- 
tude comes into excellent agreement with previous calcu- 
lations. The application of our calculations to the anal- 
ysis of the PNC measurements in Cs leads to a value of 
weak nuclear charge that is about 1.5 a smaller than the 
value predicted by the standard model. While the PNC 
amplitude found in this work is practically the same as in 
0, @| , the apparent increase in deviation from the stan- 
dard model (1 CT in when proper values of the Breit, 
radiative and neutron skin corrections are added) is due 
to smaller uncertainty. This smaller uncertainty is mostly 
due to small uncertainty of the main term which we have 
taken from Ref. [l3l without re-analysis. 

The PNC amplitude calculated in this work gives new 
constraints on physics beyond the standard model. 



CALCULATIONS 

The PNC amplitude of an electric dipole transition 
between the 6s and 7s states of cesium can be written as 
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where d — — e Ti is the electric dipole operator, Hp^c 
is the operator of a P-odd CP-even weak interaction. 
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Gf ~ 2.2225 X 10"-^^ a.u. is the Fermi constant of the 
weak interaction, and Qw is the nuclear weak charge. 

Expression ((!]) is exact if states 6s, 7s, npi/2 label 
many-electron physical states of the atom. Then 6s is 
the ground state and summation goes over excited states 
of the opposite parity and the same total angular mo- 
mentum J — 1/2. In practical calculations equation ([T]) 
is reduced to one with single-electron orbitals and single- 
electron matrix elements. It looks very similar to ([1]) 
but with a few important differences, (a) All states (6s, 
7s npi/2) are now single-electron states obtained with 
the use of the Hartree-Fock method, (b) Many-electron 
effects are reduced to redefinition of the single-electron 
orbitals and interaction Hamiltonians. For example, the 
inclusion of the core polarization effect leads to redefini- 
tion of the interaction Hamiltonian. For the weak inter- 
action we have -ffpNc = -ffpNC + <5T^nC) where SVp^c 
is the correction to the self-consistent potential of the 
atomic core due to the effect of weak interaction iJpNc- 
For the electric dipole interaction we have similar expres- 
sion d' = d + SVd. (c) Summation in ([T]) now goes over 
the complete set of single-electron states including states 
in the core. Extending summation to the core states 



corresponds to inclusion of highly excited autoionizing 
states, (d) The expression ([T|) via single-electron states 
is approximate. Its accuracy is determined by how the 
many-body effects are included. 

Equation ([1]) implies the sum-over-states approach 
which we are going to study in this paper. As mentioned 
above, high accuracy for the leading terms does not guar- 
antee high accuracy for the sum. To test the total sum 
we use an alternative approach which we have used in 
our previous PNC calculations 0, . This approach is 
based on the solving of differential equations. 

The PNC amplitude ^ can be rewritten as 



7s > 



(3) 



where the ip and Sip are single-electron orbitals and Sipa 
is the correction to the wave function ipa due to the weak 
interaction 
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It is easy to see that this correction to the wave function 
satisfies the differential equation 



(^0 - Ea)5lpa = -i?pNcV'a- 



(5) 



Here ipa is the eigenstate of the Ho Hamiltonian, which 
is in our case the relativistic Hartree-Fock (RHF) Hamil- 
tonian. The equations ([5]) have a form of the RHF equa- 
tions with the right-hand side. Solving differential equa- 
tion (O for the 6s and 7s states of cesium and using ^ 
to calculate the PNC amplitude does not require a com- 
plete set of single-electron states. It is usually numeri- 
cally more accurate than the use of the sum-over-states 
approach. In present work we use it as an independent 
test of the calculations. 

To perform the summation in ^ we use the B-spline 
basis set first presented in Ref. pjj]. We use 100 B-splines 
in each partial wave in the cavity of radius 75 as . The 
cavity radius is taken to be the same as in Ref. (lo| . 
Its value is dictated by the need to have the dominat- 
ing states be as close to physical (spectroscopic) states 
as possible. The most important intermediate states, ac- 
cording to are the 6pi/2, '^Pi/2, ^Pi/2i and 9pi/2 
states. The value i?max = 75aB is large enough for the 
9pi/2 to be physical state. The number of splines is cho- 
sen to be sufficiently large to saturate the summation. It 
turns out that saturation is achieved at approximately 
80 B-splines (in the authors used 40 B-splines of a 
different type). 

To compare the tail terms in different calculations, the 
basis sets must satisfy two conditions. First, the box size 
and number of splines must be large enough for accurate 
approximation of the all atomic states entering main term 
so that these states can be associated with real physical 
states. Second, the basis must be complete. For all basis 
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TABLE I: Partial 
10-'\{-Qw/N) a.u.] 



contributions to the -Bpnc 
for Cs in different approximations. 



[in TABLE IL Contributions of the core states npi/2 to the -EpNC 
[in 10^^^ i{ — Qw /N) a.u.] for Cs in different approximations. 



Approximation 


Core 


Main 


Tail 


Total 








Sum 


RPA" 


0.0026 


0.8705 


0.0f92 


0.8923 


DHF 


0.0247f 


-0.02645 


-0.00f74 


bo(e(2))'' 


0.00f4 


0.86f2 


0.0273 


0.8897 


RPA" 


0.0599f 


-0.0582f 


0.00f70 


bo(ae(2')'= 


O.OOfS 


0.8709 


0.0244 


0.897f 


RPA*- 


0.06043 


-0.05784 


0.00259 


bo(e'°°')'' 


O.OOfS 


0.87ff 


0.0238 


0.8967 


BO+RPA'' 


0.0723f 


-0.07049 


0.00f82 


bo(ae(°°))'= 

Ref. [fO]^ 


O.OOfS 
-0.0020 


0.8678 
0.8823 


0.0242 
0.0f95 


0.8938 
0.8998 


"RPA equations solved at Hartree-Fock frequency. 
''RPA equations solved at experimental frequency, = 


0.0844 a.u 



"Core polarization but no correlations beyond it. 

''Brueckner orbitals (BO) calculated with the second-order S. 

'^BO calculated with rescaled second-order S. 

■^BO calculated with the all-order S. 

"^BO calculated with rescaled all-order S. 

•'^DHF for the core term; coupled-cluster for the main term. 



sets which satisfy both these conditions the tail does not 
depend on basis. We believe that both our basis and that 
used in [lol| satisfy these conditions. 

We include two types of the correlations, the core po- 
larization effect and Brueckner type correlations [12| . 
The core polarization is the effect of the change in the 
self-consistent Hartree-Fock potential due to external 
field. In our case we have two types of external fields, the 
electric dipole field of the external photon and the weak 
interaction of atomic electrons with the nucleus. As we 
mentioned above, the inclusion of the core polarization 
is reduced to a redefinition of the interaction Hamiltoni- 
ans (plus small "electroweak" corrections considered in 
0, [^, It is done in the framework of the random 

phase approximation (RPA). 

The Brueckner type correlations describe the cor- 
relation interaction of the external electron with the 
atomic core, which can be reduced to redefinition of the 
single-electron orbitals, constructing the Brueckner or- 
bitals 12| . For this purpose we calculate correlation po- 



tential S |12l . |13| and construct linear combinations of 
B-splines which are eigenstates of the Hq + S Hamilto- 
nian. 



RESULTS AND DISCUSSION 

Table U shows contributions to the PNC amplitude in 
Cs obtained in different approximations within the sum- 
over-states approach. Partial contributions will be con- 
sidered in detail below. Here we start our discussion with 
some general remarks. The total PNC amplitude in the 
RPA approximation obtained with the direct summation 
(DS) is 

E^^^iBS) = 0.89235 x 10-^h{-Qw/N). 

The total PNC amplitude in the this approximation ob- 
tained via solving equations (SE) ([5]) and then using for- 



mula ([3]) is 



0.89234 X 10-^h{~Qw/N). 



Note the excellent agreement between these two num- 
bers. Since both methods of calculation share very little 
in common, it is safe to assume that numerical error is 
negligible in both cases and the number truly represents 
the PNC amplitude in the RPA approximation. This im- 
plies high quality of the basis used in present work and 
its suitability to study partial contributions. 

Table U shows significant effect of correlations on the 
tail contribution and significant difference in the val- 
ues of core and tail contributions between present work 
and There is also 1.3% difference for the main 

term. However, this difference is not important. Accu- 
rate treatment of the main term goes beyond the scope 
of the present paper; we just take its value from If 
core and tail terms also have this 1.3% relative uncer- 
tainty then it would be more than satisfactory for the 
purpose of the present work. We will see below that ac- 
tual uncertainty is higher. 

Core contribution 

Calculations in the present work are done in two dif- 
ferent ways. One uses the sum-over-states approach and 
limits the summation in ((T)) to core states. Another uses 
the weak corrections to the 6s and 7s states found by 
solving differential equations (see formula The con- 
tribution of the core states are found by imposing orthog- 
onality conditions for S'ipes and Si/jts and the states in the 
core. Both methods give exactly the same results. 

To test the calculations even further we ran our code 
for the PNC in Ra+ and found excellent agreement with 
[3] for the core contribution in RPA approximation. 

Our results are presented in Table HH The result in 
the Dirac-Hartree-Fock (DHF) approximation is in good 
agreement with those of [13 and [i] (see also Table U 
which were calculated in the DHF approximation in both 
works jl5l] . 

Note the strong cancelation between terms with Sipss 
and Sip-^s- This cancelation makes the core contribution 
very sensitive to the inclusion of the core polarization 
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effect. We include it by solving the RPA equations for 
both operators (-ffpNC and d). The equations for the 
electric dipole operator are solved a.t lu — 0.0844 a.u. 
which is the experimental energy difference between the 
6s and 7s states of Cs. 

The last line of Table [IT] presents the effect of using 
Brueckner orbitals for the core contribution. The use of 
the Brueckner orbitals in the core can be justified by the 
condition that core and valence states must be orthogonal 
to each other and using the same S operator in both 
cases is a good way to achieve this. The difference in the 
core contribution using RPA and Brueckner orbitals is 
relatively small (see Table |l| . We use this difference as 
an estimate of the uncertainty of the core contribution. 

The final difference between the present work and 
for the core contribution is 0.0038 in units of Table HIl 
This difference is mostly due to the core polarization ef- 
fect. 



Tail contribution 

The third row of Table U shows partial contributions 
to the tail component of the PNC amplitude calculated 
in different approximations. To include correlations we 
use four different sets of Brueckner orbitals obtained with 
the use of two different correlation-correction operators 
S: the second-order operator S^^^ and the all-order 
operator E(°°) 13 |. Rescaling of S is done to fit the 
energies of the lowest valence states. Rescaling usually 



improves the wave functions and therefore the matrix 
elements. 

Inclusion of the correlations leads to significant in- 
crease in the values of the tail contribution (see Table H]). 
For our final number we take our most complete calcula- 
tion, using the all-order S operator, while the spread of 
values in the various Brueckner approximations gives a 
reliable estimate of the uncertainty in our methods. 

The result of for the tail (0.0195, see Table H was 
obtained using a blend of many-body approximations in- 
cluding a simplified coupled-cluster method which only 
includes single and double excitations [15|. The result is 
close to our result in the RPA approximation but signif- 
icantly smaller than our correlated value. 



Ill: Correlated PNC amplitude in Cs [in 
Qw/N) a.u.] in different calculations. Breit, 



TABLE 

io-"j(- 

QED and neutron skin corrections are not included. 



Value 



Source and method 



0.908(9) 

0.909(9) 

0.905(9) 

0.9078(45) 

0.8998(24) 

0.9079(40) 



CP+PTSCI, Ref. [2] 
CC SD, Ref. [3] 
MBPT with fitting, Ref. [8] 
CP-hPTSCI, Ref. [9] 
CC SDvT, Ref [10] 
This work 



TABLE IV: All significant contributions to the -Epnc [in 



10~"i(-QvK/A) a.u.] for Cs. 



Contribution 



Value 



Source 



Core (n < 6) 

Main (n = 6 - 9) 

Tail (n > 9) 

Subtotal 

Breit 

QED 

Neutron skin 
Total 



0.0018(8) 
0.8823 (17) 
0.0238 (35) 
0.9079 (40) 
-0.0055 (1) 
-0.0029 (3) 
-0.0018 (5) 
0.8977 (40) 



This work 
Ref [10] 
This work 
This work 
Ref [5,, 6] 
Ref [7] ' 
Ref [5] 
This work 



ble and valence triple excitations. All numbers, apart 
from those of Ref. [lO[ are in very good agreement with 
each other. But if the result of [1011 is corrected as shown 
in Table HI it comes to very good agreement with other 
calculations as well (last line of Table IIIip . 

We summarize the results in Table IIVI We take the 
main term from Ref. [loj assuming that its value and 
uncertainty were calculated correctly. The core and tail 
contributions come from the present work. Then we add 
all other significant contributions to the PNC amplitude 
in cesium which can be found in the literature. The final 
value for the PNC amplitude is 



£^PNC = 0.8977 (40) X l{r^^i{-Qw /N) , 



(6) 



which is in excellent agreement with our previous cal- 
culations, EpNc = 0.8980(45) ^M■ The experimental 
value for the PNC amplitude is [l| 



£:pNc//3 = 1.5935 (56) mV/cm. 



(7) 



Summary 

Table [Till presents the results of the most accurate cal- 
culations of the correlated (without Breit, quantum elec- 
trodynamic (QED) and neutron skin correction) PNC 
amplitude in Cs. The abbreviation CP-I-PTSCI stands 
for correlation potential 12| combined with the pertur- 
bation theory in screened Coulomb interaction, CC SD 
stands for coupled cluster with single and double excita- 
tions, CC SDvT means coupled cluster with single, dou- 



The most accurate value for the vector transition proba- 
bility /3 comes from the analysis of the Bennett and 
Wieman measurements 1171 



(3 = 26.957 (51) a|. 
Comparing ([6|), jZl) and ([8|) leads to 



(8) 



Q^("3CS) = ^72.58 (29)expt (32)theory • (9) 

This value is in a reasonable agreement with the pre- 



diction of the standard model, = -73.23(2) [18 
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If we add theoretical and experimental errors in ([9]) in 
quadrature, the Cs PNC result deviates from the stan- 
dard model value by 1.5 a: 



Qm' - g|f = (JQw- = 0.65 (43). 



(10) 



For small deviations from the Standard Model values we 
may relate this to the deviation in sin'^ 9\y using the sim- 
ple relationship 5Qw ~ —4:Z 6{sin^ 0\y) which gives 



sm 



0.2356(20). 



(11) 



value 



This is a|so 1.5 a off the standard model 
0.2386 (1) [l^ at near zero momentum transfer. 

The new physics originated through vacuum polariza- 
tion to gauge boson propagators is described by weak 
isospin conserving S and isospin breaking T parame- 
ters [il 



Qw - Qlf = -0.800 5- 0.007 T. 



(12) 



At the 1(7 level leads to 5 = -0.81 (54). 

Finally, a positive AQw could also be indicative of an 
extra Z boson in the weak interaction 1201 



■0A{2N + Z){Mw/Mzx) 



(13) 



Using ([To]) leads to Mz^ > 710 GeV/c^. 

The authors are grateful to A. Derevianko for useful 
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tralian Research Council. 
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